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LONG-TERM  GOALS 

Sea  ice  reflectance  or  albedo,  a  key  parameter  in  climate  modeling,  is  primarily  determined  by  melt 
pond  and  ice  floe  configurations.  Ice-albedo  feedback  has  played  a  major  role  in  the  recent  declines  of 
the  summer  Arctic  sea  ice  pack.  However,  understanding  the  evolution  of  melt  ponds  and  sea  ice 
albedo  remains  a  significant  challenge  to  improving  climate  models.  Our  research  is  focused  on 
obtaining  extensive  imagery  of  melt  pond  evolution,  and  developing  mathematical  models  of  the 
melting  process  that  can  help  us  better  understand  the  role  of  sea  ice  in  the  climate  system,  and 
represent  sea  ice  more  rigorously  in  climate  models. 

OBJECTIVES 

Viewed  from  high  above,  the  melting  sea  ice  surface  can  be  thought  of  as  a  two  phase  composite  of  ice 
and  melt  water.  The  boundaries  between  the  two  phases  evolve  with  increasing  complexity  and  a  rapid 
onset  of  large  scale  connectivity,  or  percolation  of  the  melt  phase.  We  plan  to  document  this 
phenomenon  with  photographic  imagery  and  to  develop  percolation  and  other  models  to  quantitatively 
describe  the  evolution. 

APPROACH 

Key  Particpants: 

Ken  Golden  (Professor  of  Mathematics,  U.  of  Utah) 

Don  Perovich  (Research  Geophysicist,  ERDC-CRREL) 

Tolga  Tasdizen  (Associate  Professor,  Department  of  Electrical  and  Computer  Engineering,  Scientific 
Computing  and  Imaging  Institute,  U.  of  Utah) 

Court  Strong  (Assistant  Professor,  Department  of  Atmospheric  Sciences,  U.  of  Utah) 
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We  describe  here  various  facets  of  our  approach  to  characterizing  and  modeling  Arctic  melt  pond 
evolution. 

Acquisition  of  melt  pond  imagery  (Don  Perovich,  Chris  Polashenski,  and  Alexandra  Arntsen):  We 
are  developing  a  library  of  melt  pond  imagery  that  will  be  analyzed.  This  library  includes  aerial 
photographs  giving  complete  time  series  of  melt  pond  evolution  during  the  1997-1998  Surface  Heat 
Budget  of  the  Arctic  Ocean  field  campaign  plus  a  trans-Arctic  set  of  photographs  from  August- 
September  2005  showing  spatial  variability  of  ponds.  In  addition,  we  are  archiving  satellite  images 
from  Quikbird  and  National  Technical  Means  satellites. 

Mapping  melt  pond  images  onto  discrete  networks  (Tolga  Tasdizen,  Meenakshi  Barjatia,  Boya 
Song,  Ken  Golden):  Given  a  photograph  of  a  complex  melt  pond,  we  are  developing  algorithmic 
methods  for  mapping  the  melt  pond  onto  a  resistor  network  so  that  we  can  study  the  evolution  of 
connectivity,  look  for  a  percolation  threshold,  and  compute  horizontal  flow  characteristics. 


Modeling  projects 

We  are  developing  mathematical  characterizations  of  melt  pond  evolution  using  models  of  composite 
materials  and  statistical  physics.  In  earlier  work  done  just  prior  to  this  grant,  we  discovered  that,  as 
melt  ponds  grow,  coalesce,  and  become  more  complex,  they  exhibit  a  transition  in  fractal  dimension 
from  about  1  to  2,  near  a  critical  length  scale  of  100  square  meters.  We  have  been  pursuing  a  number 
of  avenues  to  model  the  evolution  of  melt  pond  configurations  and  better  understand  the  transition  in 
complexity  that  real  ponds  display. 

Continuum  percolation  model  (Brady  Bowen,  Court  Strong,  and  Ken  Golden):  We  have  developed  a 
continuum  percolation  model  of  melt  pond  evolution.  The  sets  created  by  the  intersection  of  a  plane 
(water  level)  with  a  surface  generated  by  a  random  Fourier  series  (representing  the  snow  and  ice 
topography)  look  very  similar  to  melt  ponds,  as  the  plane  is  moved  up  and  down. 

Ising  model  (Ivan  Sudakov,  Yi-Ping  Ma,  and  Ken  Golden):  The  Ising  model  has  been  perhaps  the 
most  widely  studied  model  of  a  phase  transition  in  statistical  physics.  The  model  describes  a  lattice  of 
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magnetic  spins  that  interact  with  an  applied  magnetic  field  as  well  as  with  each  other.  We  are  working 
to  adapt  this  model  to  describe  melt  pond  evolution.  The  key  magnetic  states,  spin  up  and  spin  down, 
are  replaced  by  melt  water  and  ice. 

Network  model  of  horizontal  fluid  permeability  (Boya  Song  and  Ken  Golden):  Drainage  of  melt 
ponds  can  occur  not  only  through  porous  sea  ice  below  the  pond,  but  at  a  distance  though  seal  holes, 
cracks,  etc.  Horizontal  flow  through  fluid  pathways  facilitates  this  type  of  drainage.  We  have  been 
working  on  creating  resistor  network  models  of  horizontal  fluid  flow,  and  computing  the  horizontal 
permeability  of  melt  pond  configurations. 

Spectral  measures  for  melt  pond  configurations  (Ben  Murphy  and  Ken  Golden):  The  analytic 
continuation  method  provides  rigorous  integral  representations  for  the  effective  transport  coefficients 
of  composite  media.  The  key  mathematical  object  in  an  integral  representation  is  a  so-called  “spectral 
measure”  which  depends  only  on  the  composite  geometry.  In  the  previous  sub-project  the  importance 
of  horizontal  fluid  transport  is  described,  and  computation  of  these  spectral  measures  provides  a 
general  technique  for  computing  horizontal  flow  properties,  as  well  as  other  parameters  of  interest. 

Voter  model  and  the  fractal  geometry  of  melt  ponds  (Rebecca  Nickerson  and  Ken  Golden): 
Interacting  particle  systems  have  been  used  to  model  the  evolution  of  a  broad  range  of  systems.  We 
have  developed  a  simple  voter-type  model  to  study  how  melt  ponds  fonn.  The  key  step  is  to  mimic  the 
assumption  in  the  voter  model  that  a  voter  may  be  more  likely  to  vote  for  a  candidate  if  his  or  her 
neighbors  have.  For  melt  ponds,  a  square  is  more  likely  to  melt  if  more  of  its  neighbors  have  melted. 

Cheeger  constants  for  melt  ponds  (Kasey  Leavitt  and  Ken  Golden):  We  have  been  developing 
methods  of  algorithmically  mapping  images  of  melt  ponds  onto  random  planar  graphs.  To  measure  key 
geometrical  features  of  the  ponds,  we  have  been  investigating  how  to  compute  so-called  Cheeger 
constants  for  graphs  representing  melt  pond  configurations.  This  parameter  is  a  measure  of  whether  or 
not  the  graph  has  a  “bottleneck,”  which  is  a  critical  geometrical  feature  controlling  horizontal  flow  of 
melt  water.  The  notion  of  the  Cheeger  constant  should  also  be  useful  for  analysis  of  the  brine 
microstructure  of  sea  ice. 

Melt  Ponds  and  Bifurcations  in  Climate  Models  (Ivan  Sudakov  and  Ken  Golden):  As  melt  ponds 
grow  and  coalesce,  their  fractal  dimension  changes,  which  means  that  their  areas  don’t  grow  as  rapidly 
as  their  perimeters.  Does  this  impact  albedo  and  the  dynamics  of  conceptual  climate  models?  We  are 
investigating  this  possibility. 

Nonlinear  PDE  model  for  melt  pond  water  depth  (Court  Strong  and  Ken  Golden):  In  order  to 
provide  a  physically  based  mathematical  framework  to  complement  our  statistical  physics  analyses  of 
melt  pond  fractal  dimension,  we  have  been  working  with  a  partial  differential  equation  governing  melt 
pond  water  depth  h  (Luthje  et  al.,  2006). 

Numerical  model  for  fluid  flow  through  porous  sea  ice  (Kyle  Steffen,  Yekaterina  Epshteyn,  and 
Ken  Golden):  Melt  pond  evolution  is  intimately  related  to  how  easy  it  is  for  fluid  to  flow  or  drain 
through  the  porous  brine  microstructure  of  the  sea  ice  supporting  the  ponds.  We  have  begun 
developing  numerical  models  of  this  important  process,  and  modifying  existing  random  pipe  models  to 
investigate  the  effect  of  the  presence  of  algal  EPS  on  fluid  penneability. 
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Field  Project 

Chris  Polashenski  and  Ken  Golden  participated  in  the  SUB  ICE  (Study  of  Under  ice  Blooms  In  the 
Chuckchi  Ecosystem)  Expedition  aboard  the  USCGC  Healy  for  about  six  weeks  during  May  and  June 
of  2014  (Don  Perovich  was  the  leader  of  the  ice  team  but  was  unable  to  travel  due  to  a  medical  issue, 
Polashenski  took  over  the  leadership  role).  We  looked  at  ice  albedo  and  melt  pond  fonnation. 


WORK  COMPLETED 

Image  Acquisition.  We  are  building  a  library  of  aerial  and  satellite  imagery.  So  far  we  have  compiled 
a  few  thousand  aerial  photographs  of  ponded  sea  ice  plus  approximately  50  satellite  images. 

Fractal  dimension  computation.  In  discovering  that  there  is  a  transition  in  the  fractal  dimension  of 
melt  ponds,  we  needed  to  come  up  with  a  way  of  computing  the  fractal  dimension  from  area-perimeter 
data  on  melt  ponds.  We  developed  a  detenninistic  algorithm  to  do  this.  However,  a  statistical  approach 
that  creates  a  “best  fit”  fractal  dimension  curve  for  the  given  data  set  eluded  us.  Court  Strong  found 
such  a  method.  It  creates  an  optimal  fit  of  a  hyperbolic  tangent  model  for  the  fractal  dimension  as  a 
function  of  log  A,  where  A  is  the  area  of  the  melt  ponds.  This  algorithm  now  makes  possible  the  fractal 
analysis  of  melt  pond  configurations  produced  by  mathematical  models,  as  well  as  raw  melt  pond  data 
that  we  will  be  obtaining  in  the  future. 

Melt  pond  evolution  models.  The  pond  configurations  in  the  continuum  percolation  model  are  for  a 
surface  produced  by  a  Fourier  series  with  random  amplitudes  and  phases.  These  configurations  look 
quite  similar  to  actual  melt  ponds.  Similarly,  the  voter  model  produces  configurations  that  look  like 
melt  ponds,  as  does  the  Ising  model.  The  challenge  now  is  to  better  understand  the  basic  physics  of 
melt  ponds  and  how  the  parameters  of  these  models  are  related  to  melt  pond  characteristics.  We  have 
done  this  for  an  Ising  model  of  melt  pond  evolution,  having  now  incorporated  ice-albedo  feedback  and 
the  basic  thermodynamics  of  melting  ice. 

Network  model  of  horizontal  fluid  permeability.  We  have  worked  through  examples  where  a 
complex  melt  ponds  are  mapped  onto  a  resistor  network.  This  network  contains  not  only  connectivity 
information  about  the  different  components  of  the  pond,  but  the  widths  of  narrow  channels  through 
which  melt  water  must  flow,  for  example,  to  be  drained  through  a  nearby  seal  hole.  The  penneability 
calculation  depends  on  computing  the  effective  conductivity  of  the  network.  We  have  implemented  an 
exact  fonnula  for  the  conductivity  of  a  random  graph  -  based  on  an  enumeration  of  trees  of  the  graph 
via  the  determinants  of  adjacency  matrices  -  in  order  to  compute  the  horizontal  fluid  permeability. 
Moreover,  we  have  developed  an  algorithmic  approach  to  mapping  melt  pond  images  onto  these 
networks. 

Spectral  Measures  for  Melt  Pond  Configurations.  We  have  computed  examples  of  spectral 
measures  for  discretizations  of  actual  melt  pond  images,  and  have  made  significant  advances  in  the 
underlying  functional  and  numerical  analysis  needed  for  these  computations. 

Nonlinear  PDE  model  for  melt  pond  water  depth.  We  developed  code  to  numerically  integrate  the 
melt  pond  water  depth  PDE  on  grids  with  various  horizontal  resolutions,  and  initialized  the  model  with 
topography  adapted  from  LiDAR  observations  (e.g.,  Polashenski  et  ah,  2010).  This  framework 
developed  melt  ponds  with  realistic  fractal  properties,  and  we  are  investigating  how  these  properties 
depend  on  the  specification  of  the  initial  topography. 
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Image  analysis  of  melt  ponds.  Part  of  our  effort  is  focused  on  exploring  the  evolution  of  melt  pond 
fraction,  number,  size,  and  shape  through  the  use  of  National  Technical  Means  Imagery.  This  is  a 
Lagrangian  data  set  of  imagery  that  tracks  individual  floes.  An  image  processing  technique  has  been 
developed  to  analyze  melt  pond  formation  and  evolution  during  the  summer  melt  season.  Statistics  can 
be  extracted  for  an  ensemble  of  ponds  developing  on  an  individual  floe  from  the  formation  of  discrete 
melt  ponds  to  a  network  of  connected  melt  features.  From  the  processed  imagery,  pond  area  fraction, 
pond  count,  and  average  pond  size  can  be  calculated  at  several  times  throughout  the  season. 

This  method  for  image  processing  begins  with  an  object  based  classification  that  includes  a 
segmentation  of  the  original  image  into  features  of  similar  spectral,  spatial,  and  texture  characteristics 
using  the  ENVI  Feature  Extraction  workflow.  A  series  of  variables  are  extracted  for  each  feature  and 
construct  a  shapefile  of  polygons  each  with  21  attributes.  A  histogram  analysis  and  stretch  are 
performed  to  create  additional  attributes  indicating  pixel  value  relative  to  the  entire  image.  A  training 
set  is  then  created  from  four  different  images  varying  in  time  of  season  and  lighting  conditions  for 
classes:  ice,  open  water,  submerged  ice,  and  ridges.  A  random  forest  classification  is  perfonned  in 
Matlab  to  classify  each  segmented  image  using  the  training  set.  The  result  is  a  shapefile  of  classified 
polygons.  The  class  called  submerged  ice  is  separated  into  ponds  and  submerged  ice  along  floe  edges 
using  the  spatial  proximity  to  open  water  ARCGIS  neighbor  analysis. 


Figure  1.  Images  from  June  5  and  June  18,  2013.  The  left  column  shows  the  raw  image  and  the  right 
column  the  processed  image. 


5 


RESULTS 


Below  we  present  abstracts  from  papers  that  contain  results  from  some  of  the  projects  described  above. 

Arctic  melt  ponds  and  bifurcations  in  the  climate  system 

I.  Sudakov,  S.  A.  Vakulenko  and  K.  M.  Golden 

Abstract:  Understanding  how  sea  ice  melts  is  critical  to  climate  projections.  In  the  Arctic,  melt  ponds 
that  develop  on  the  surface  of  sea  ice  floes  during  the  late  spring  and  summer  largely  determine  their 
albedo  -  a  key  parameter  in  climate  modeling.  Here  we  explore  the  possibility  of  a  conceptual  sea  ice 
climate  model  passing  through  a  bifurcation  point  -  an  irreversible  critical  threshold  as  the  system 
warms,  by  incorporating  geometric  infonnation  about  melt  pond  evolution.  This  study  is  based  on  a 
bifurcation  analysis  of  the  energy  balance  climate  model  with  ice-albedo  feedback  as  the  key 
mechanism  driving  the  system  to  bifurcation  points. 


Ising  model  for  melt  ponds  on  Arctic  sea  ice 

Y.  Ma,  I.  Sudakov,  and  K.  M.  Golden 

Abstract:  The  albedo  of  melting  Arctic  sea  ice,  a  key  parameter  in  climate  modeling,  is  detennined  by 
pools  of  water  on  the  ice  surface.  Recent  observations  show  an  onset  of  pond  complexity  at  a  critical 
area  of  about  100  square  meters,  attended  by  a  transition  in  pond  fractal  dimension.  To  explain  this 
behavior  and  provide  a  statistical  physics  approach  to  sea  ice  modeling,  we  introduce  a  two 
dimensional  Ising  model  for  pond  evolution  which  incorporates  ice-albedo  feedback  and  the 
underlying  thermodynamics.  The  binary  magnetic  spin  variables  in  the  Ising  model  correspond  to  the 
presence  of  melt  water  or  ice  on  the  sea  ice  surface.  The  model  exhibits  a  second-order  phase  transition 
from  isolated  to  clustered  melt  ponds,  with  the  evolution  of  pond  complexity  consistent  with  the 
observations. 


Network  modeling  of  Arctic  melt  ponds 

M.  Barjatia,  T.  Tasdizen,  B.  Song,  and  K.  M.  Golden 

Abstract:  The  recent  precipitous  losses  of  summer  Arctic  sea  ice  have  outpaced  the  projections  of 
most  climate  models.  Efforts  to  improve  these  models  have  focused  in  part  on  a  more  accurate 
accounting  of  sea  ice  albedo  or  reflectance.  In  late  spring  and  summer,  the  albedo  of  the  ice  pack  is 
detennined  primarily  by  melt  ponds  that  form  on  the  sea  ice  surface.  The  transition  of  pond 
configurations  from  isolated  structures  to  interconnected  networks  is  critical  in  allowing  the  lateral 
flow  of  melt  water  toward  drainage  features  such  as  large  brine  channels,  fractures,  and  seal  holes, 
which  can  significantly  alter  the  albedo.  Moreover,  pond  connectivity  can  also  influence  their 
effectiveness  in  breaking  up  an  ice  floe  as  the  melt  season  progresses.  Here  we  develop  algorithmic 
techniques  for  mapping  photographic  images  of  melt  ponds  onto  discrete  conductance  networks  which 
represent  the  geometry  of  pond  configurations  and  approximate  the  ease  of  lateral  flow.  We  implement 
an  image  processing  algorithm  with  mathematical  morphology  operations  to  produce  a  conductance 
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matrix  representation  of  the  melt  ponds.  Basic  clustering  and  edge  elimination  using  graph  theory  are 
then  used  to  reduce  the  conductance  matrix  to  include  only  direct  connections.  The  results  for  images 
taken  during  different  times  of  the  year  are  visually  inspected  and  the  number  of  mislabels  is  used  to 
evaluate  perfonnance. 


Fractal  properties  of  a  continuum  percolation  model  for  the  evolution  of  melt  ponds  on  Arctic 
sea  ice 

B.  Bowen,  C.  Strong,  and  K.  M.  Golden 

Abstract:  Melt  ponds  on  Arctic  sea  ice  are  known  to  undergo  a  dramatic  transition  in  fractal 
dimension  as  they  grow  and  connect  during  the  melt  season.  Here,  we  employ  a  continuum  percolation 
model  where  the  pond  boundaries  are  the  intersections  of  a  horizontal  plane  with  a  random  surface 
representing  ice  topography.  We  define  a  class  of  stochastic  surfaces  that  replicate  this  fractal 
dimension  shift,  and  we  explore  the  fractal  properties  of  the  associated  ponds  as  a  function  of  their 
autocorrelation  and  anisotropy.  We  use  two  models  to  investigate  the  evolution  of  melt  ponds  -  the 
level  set  percolation  model  where  the  melt  ponds  are  the  volume  between  the  rising  plane  and  the 
surface,  and  a  partial  differential  equation  model  based  on  water  mass  conservation  and  accelerated 
melt  rates  due  to  ice-albedo  feedback  where  water  is  present.  Both  models  exhibit  a  transition  in  the 
fractal  dimension,  but  we  find  that  the  partial  differential  equation  model  shifts  to  a  larger  fractal 
dimension.  The  water-enhanced  melt  rates  deepen  rather  than  widen  the  ponds,  allowing  them  to  grow 
and  connect  which  increases  perimeter  without  substantial  expansion  of  their  surface  area. 


Random  matrix  universality  for  connectedness  transitions  in  composite  materials 

N.  B.  Murphy  and  K.  M.  Golden 

Abstract:  Random  matrix  universality  in  large  complex  systems,  where  the  eigenvalues  of  random 
matrix  models  display  universal  statistical  behavior,  arises  from  a  fine  balance  between  chaos  and 
order.  Universality  has  turned  up  in  so  many  seemingly  unrelated  problems  that  it  can  be  thought  of  as 
a  fundamental  law,  not  unlike  the  central  limit  theorem.  A  generalization  of  classical  universality  has 
been  successful  in  describing  some  chaotic  quantum  systems  and  disordered  mesoscopic  conductors 
which  undergo  transitional  behavior  as  system  parameters  vary.  We  find  that  this  generalized 
universality  also  governs  transitional  behavior  in  the  effective  transport  properties  of  macroscopic 
composite  materials.  In  particular,  transport  in  composites  displays  critical  behavior  near  transitions  in 
the  connectedness  or  percolation  properties  of  a  particular  phase.  The  analytic  continuation  method 
provides  rigorous  integral  representations  for  the  bulk  transport  coefficients  of  such  media.  At  the  heart 
of  the  method  is  a  family  of  random  matrices,  new  to  random  matrix  theory,  which  depend  only  on  the 
composite  geometry.  Here  we  demonstrate  that  an  onset  of  connectedness  gives  rise  to  striking 
transitional  behavior  in  the  long  and  short  range  correlations  of  the  eigenvalues  of  the  random 
matrices.  Their  statistics  shift  from  weakly  to  highly  correlated,  or  Poissonian  toward  Wigner-Dyson, 
as  a  function  of  system  connectedness.  This  provides  a  mechanism  for  explaining  the  collapse  of 
spectral  gaps  which  control  critical  behavior  of  transport.  The  spectral  transition  is  explored  for 
random  resistor  networks,  human  bone,  and  some  sea  ice  structures  with  a  range  of  scales,  whose 
effective  properties  are  important  in  climate  modeling. 
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Statement  of  significance:  Since  the  1950s  during  studies  of  atomic  spectra,  the  striking  phenomenon 
known  as  random  matrix  universality  has  emerged  in  the  signal  dynamics  of  the  Internet,  the  departure 
times  for  bus  systems,  and  even  in  the  zeros  of  the  famous  Riemann  zeta  function.  The  eigenvalues  or 
energy  levels  of  random  matrix  models  of  such  systems  exhibit  universal  statistical  behavior.  A 
generalized  universality  has  accurately  described  transitions  in  quantum  and  mesoscopic  systems.  We 
discover  here  that  this  generalized  universality  also  describes  transitions  in  the  connectedness  or 
percolation  properties  of  composite  materials  with  microstructural  scales  ranging  from  microns  to 
kilometers.  The  associated  transport  properties  of  these  composites  arise  in  many  engineering  and 
technological  applications,  and  play  a  key  role  in  Earth's  climate  system. 


Spectral  measure  computations  for  composite  media 

N.  B.  Murphy,  C.  Hohenegger,  E.  Cherkaev  and  K.  M.  Golden 

Abstract:  The  analytic  continuation  method  of  homogenization  theory  provides  Stieltjes  integral 
representations  for  the  effective  parameters  of  composite  media.  These  representations  involve  the 
spectral  measures  of  self-adjoint  random  operators  which  depend  only  on  the  composite  geometry.  On 
finite  bond  lattices,  these  random  operators  are  represented  by  random  matrices  and  the  spectral 
measures  are  given  explicitly  in  tenns  of  their  eigenvalues  and  eigenvectors.  Here  we  provide  the 
mathematical  foundation  for  rigorous  computation  of  spectral  measures  for  such  composite  media,  and 
develop  a  numerically  efficient  projection  method  to  enable  such  computations.  This  is  accomplished 
by  providing  a  novel  fonnulation  of  the  analytic  continuation  method  which  is  equivalent  to  the 
original  formulation  and  holds  for  both  the  finite  lattice  setting  and  the  infinite  settings.  We  also 
introduce  a  family  of  random  bond  lattices  and  directly  compute  the  associated  spectral  measures  and 
effective  parameters.  The  computed  spectral  measures  are  in  excellent  agreement  with  known 
theoretical  results.  The  behavior  of  the  associated  effective  parameters  is  consistent  with  the 
symmetries  and  theoretical  predictions  of  models,  and  the  computed  values  fall  within  rigorous 
bounds.  Some  previous  calculations  of  spectral  measures  have  relied  on  finding  the  boundary  values  of 
the  imaginary  part  of  the  effective  parameter  in  the  complex  plane.  Our  method  instead  relies  on  direct 
computation  of  the  eigenvalues  and  eigenvectors  which  enables,  for  example,  statistical  analysis  of  the 
spectral  data. 


Freshwater  -  the  key  to  melt  pond  formation  atop  first  year  sea  ice 

C.  Polashenski,  K.  M.  Golden,  E.  Skyllingstad,  and  D.  Perovich 

Abstract:  Melt  pond  formation  atop  Arctic  sea  ice  is  a  primary  control  of  shortwave  energy  balance 
and  light  availability  for  photosynthesis  in  the  upper  Arctic  Ocean.  The  initial  fonnation  process  of 
melt  ponds  on  first  year  ice  typically  requires  that  melt  water  be  retained  on  the  surface  of  ice  several 
to  tens  of  centimeters  above  sea  level  for  several  days.  Albedo  feedbacks  during  this  time  period  create 
below-sea-level  depressions  which  remain  ponds  later  in  summer.  Both  theory  and  observations, 
however,  show  that  sea  ice  is  so  highly  porous  and  permeable  prior  to  the  formation  of  melt  ponds  that 
retention  of  water  tens  of  cm  above  hydraulic  equilibrium  for  multiple  days  should  not  be  possible. 
Here  we  present  results  of  percolation  test  experiments  that  identify  the  mechanism  allowing  above-sea 
level  melt  pond  fonnation.  The  infiltration  of  fresh  water  into  pore  structure  of  the  ice  is  responsible 
for  plugging  the  pores  with  fresh  ice,  sealing  the  ice  against  further  water  percolation,  and  allowing 
water  to  pool  above  freeboard.  Fresh  meltwater  availability  and  desalination  processes,  therefore, 
likely  exert  considerable  influence  over  the  fonnation  of  melt  ponds.  The  findings  demonstrate  another 
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mechanism  through  which  changes  in  snowfall  on  sea  ice,  already  being  observed,  are  likely  to  alter 
ice  mass  balance  and  highlight  the  importance  of  efforts  to  improve  treatment  of  ice  salinity  in  models. 


Object-based  detection  of  Arctic  sea  ice  and  melt  ponds  using  high  spatial  resolution  aerial 
photographs 

X.  Miao,  H.  Xieb,  S.  F.  Ackley  and  D.  K.  Perovich 

Abstract:  High  resolution  aerial  photographs  used  to  detect  and  classify  sea  ice  features,  especially 
melt  ponds,  can  provide  extracted  physical  parameters  to  refine,  validate,  and  improve  climate  models. 
However,  manually  delineating  sea  ice  and  melt  ponds  is  time-consuming  and  labor-intensive.  In  this 
study,  an  object-based  classification  scheme  is  used  to  extract  sea  ice  and  melt  ponds  efficiently  from 
163  selected  aerial  photographs  taken  during  the  Chinese  National  Arctic  Research  Expedition 
(CHINARE  2010)  to  the  western  Arctic.  The  algorithm  includes  three  major  steps:  (1)  the  image 
segmentation  groups  the  neighboring  pixels  into  objects  according  to  the  similarity  of  spectral  and 
textural  information;  (2)  the  random  forest  ensemble  classifier  distinguishes  four  general  classes: 
water,  general  submerged  ice  (GSI),  shadow,  and  ice/snow;  and  (3)  the  polygon  neighbor  analysis 
further  separates  melt  ponds  and  submerged  ice  from  GSI  according  to  their  spatial  relationships.  The 
overall  classification  accuracy  for  the  four  general  classes  is  95.5%  based  on  178  ground  reference 
objects.  Furthermore,  the  producer’s  accuracy  of  90.8%  and  user’s  accuracy  of  91.8%  is  achieved  for 
melt  pond  detection  through  98  independent  reference  objects.  For  the  163  photos  in  the  Arctic 
marginal  ice  zone  examined,  a  total  of  19,438  melt  ponds  larger  than  1  m  are  detected,  with  a  pond 
density  of  867.2  km-1.  The  average  pond  size  is  32.6  m2,  and  the  average  melt  pond  fraction  is  5.63%. 
This  method  can  be  applied  to  massive  high  spatial  resolution  Arctic  sea  ice  photographs  for  deriving 
detailed  sea  ice  and  melt  pond  distributions  over  wider  regions,  and  extracting  sea  ice  physical 
parameters  and  their  corresponding  changes  between  years. 


IMPACT/APPLICATIONS 

The  simplicity  of  some  of  our  melt  pond  models  holds  out  the  possibility  of  efficiently  accounting  for 
melt  ponds  and  sea  ice  albedo  in  climate  models.  Our  Ising  model  represents  the  first  application  of  a 
sophisticated  statistical  physics  model  to  studying  the  evolution  of  melt  ponds.  Our  field  work  answers 
the  fundamental  puzzle  about  how  melt  ponds  can  form  and  persist  atop  highly  permeable  sea  ice. 


RELATED  PROJECTS 

NSF  DMS  Math  Climate  Research  Network  Grant  (http://www.mathclimate.org/).  This  grant  is 
partially  funding  our  melt  pond  work  and  research  on  related  sea  ice  properties. 
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